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Abstract 

We study the thermal leptogenesis in the Eq x 17(1) a SUSY GUT 
model in which realistic masses and mixings of quarks and leptons can 
be realized. We show that the sufficient baryon number can be pro¬ 
duced by the leptogenesis in the model, in which the mass parameter 
of the lightest right-handed neutrino is predicted to be smaller than 
10^ GeV. The essential point is that the mass of the lightest right- 
handed neutrino can be enhanced in the model because it has a lot of 
mass terms whose mass parameters are predicted to be the same or¬ 
der of magnitude which is smaller than 10® GeV. We show that 0(10) 
enhancement for the lightest right-handed neutrino mass is sufficient 
for the observed baryon asymmetry. Note that such mass enhance¬ 
ments do not change the predictions of neutrino masses and mixings 
at the low energy scale in the Eq model which has six right-handed 
neutrinos. In the calculation, we include the effects of supersymme¬ 
try and ffavor in final states of the right-handed neutrino decay. We 
show that the effect of supersymmetry is quite important even in the 
strong washout regime when the effect of ffavor is included. This is 
because the washout effects on the asymmetries both of the muon and 
the electron become weaker than that of the tau asymmetry. 
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1 Introduction 


Supersymmetric (SUSY) grand unified theory (GUT)[T] is one of the most 
promising candidates as the extended model of the standard model (SM). 
This is because the SUSY GUT realizes two kinds of unihcations, unihcation 
of the gauge interactions and unification of the matters in the SM and for both 
unihcations, there are supports from experiments. Three gauge couplings in 
the SM meets at a scale, which is called the GUT scale Aq ~ 2 X 10^® GeV. 
Moreover, the various hierarchies of quark and lepton masses and mixings can 
be naturally understood in SU (5) unihcation if we assume that the 10 helds 
of SU{5) induce stronger hierarchy in Yukawa couplings than the 5 helds 
of SU{5). One of the most important advantages of the Eq unihcation|2] is 
that the above assumption can be naturally derived[3]. As the result of this 
important feature of the Eq unihcation, we can build an Eq GUT in which 
all three generation of quarks and leptons can be unihed into a single multi- 
plet(or two multiplets) by introducing family symmetry SU{3)p{ot SU{2)p) 
and the realistic quark and lepton masses and mixings can be realized after 
breaking the family and GUT symmetries [1]. 

However, it is well-known that SUSY GUTs are suhering from the doublet- 
triplet splitting problem [5]. The doublet Higgs must have the weak scale 
mass to obtain the weak scale, while the triplet (colored) Higgs which be¬ 
longs to the same multiplet as the doublet Higgs in the GUT must have 
the GUT scale mass to stabilize the nucleon. Fortunately, if the anoma¬ 
lous U{1)a gauge symmetry[6] is introduced, the problem can be solved un¬ 
der a natural assumption that all the interactions are introduced with 0(1) 
coefficients [3 El E]. Because of this natural assumption, the coefficients of 
the terms and the vacuum expectation values (VEVs) of the GUT Higgs can 
be determined only by the symmetry of the theory. The coefficients of the 
interaction XYZ are determined [10 1 ITT] except the 0(1)coefficients by the 
total anomalous U (1)^ charge x + y + z as 

X^+y+^XYZ {x-\-yz > Q) . . 

0 {x -\-y + z <Q), ^ 

where x, y, and x are the U{1)a charges of the helds Y, Y, and Y, respec¬ 
tively. Throughout this paper, we denote all the helds with uppercase letters 
and their anomalous U{1)a charges with the corresponding lowercase letters 
if there is no special comment. Here A is the ratio of the Fayet-Illiopoulos 
parameter ^ to the cutoh A, and in this paper we take A ~ 0.22 as a typical 
value. Under the natural assumption, we can obtain the realistic Yukawa 
couplings in Eq GUT[3] (or in 50(10) GUT[7] which has similar structure as 
Eq gut). The VEVs of the operators O is also determined[3] by their total 
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anomalous U{1)a charges o as 


(O) = / *^ ^ (2) 

\ A-° (o < 0) • ^ ^ 

In this paper, we often use a unit in which the cutoff A is taken to be 1. 
Because of the natural assumption, all the mass spectrum of superheavy 
particles and the VEVs of GUT Higgs are determined only by the symmetry 
of the theory. Therefore, we can calculate the running gauge couplings once 
we £x the symmetry of the theory. Interestingly, this natural scenario gives a 
novel explanation [H] for the experimental support for the unihcation of three 
gauge interactions in the SM. The new explanation requires that the cutoff 
scale must be taken to be around the usual GUT scale Ac [9]. 

If this natural Eq GUT describes our world, it must be consistent with 
the cosmology. The dark matter can be the lightest supersymmetric parti¬ 
cle. In this paper, we discuss the leptogenesis|12j in this scenario. One of 
the important things in unihcation for the leptogenesis is that the funda¬ 
mental representation 27, which is decomposed in the Eq D S'O(IO) x U{l)v' 
notation (and in the [50(10) D SU{5) x U{l)v] notation) as 

27 = 16i[10i -|- 5_3 -|- I 5 ] -|- 10_2[5_2 -I- -I- (3) 

includes two singlets 5(1') and A^r( 1) under the SM gauge group, which 
can be the right-handed (RH) neutrinos. If we introduce three 27 for three 
generation quarks and leptons, we have six RH neutrinos. Basically, since 
the masses and Yukawa couplings of the RH neutrinos are determined by the 
symmetry, we can examine whether the leptogenesis works well or not in this 
scenario. Naively, the leptogenesis in this scenario does not work because 
the lightest RH neutrino becomes lighter than 10® GeV, i.e., this scenario 
looks not to satisfy the Ibarra’s upper bound [13] for the lightest RH neu¬ 
trino which is 10 ®“® GeV. Actually, in a typical model, the (effective) U{1)a 
charges of Si and {i = 1,2,3) are hxed as (si,S 2 ,S 3 ) = (6.5, 5.5, 3.5) 
and (h^i,h ^2 5 = (6,5,3), and therefore, the mass of the lightest RH 
neutrino 5i becomes Ms^ Ai®A ~ 5.7 X 10^ GeV[3]. Yukawa couplings 
are also easily estimated because the sum of the (effective) U{1)a charges of 
the up-type Higgs Hu and doublet-leptons k become {hu + h, hu + h, hu + 
h) = (0,—0.5,—1). The Yukawa couplings among /*, 5i and Hu become 
(A®■^A^A5■5). Then, we can estimate two important parameters for the lep¬ 
togenesis as 

K = Td/H ~ 40 

g ^ r{Si^i+Hu)-riSi^I+Hl) ^ 5 ^ ]^q-9 (4) 
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where T£> and H are the decay width of Si and the Hubble parameter at 
T = M 5 J, respectively. (In this paper we denote the lepton doublet helds 
with lowercase letter I in order to avoid the confusion with lepton asymmetry 
L in the following discussions. ) Since the sufficient production of Baryon 
number requires ~ 1 and e ~ 10“^, this K is too large, and the e is too 
small. The produced lepton number is calculated as 

Yl = — ^ 10-13^ (5) 

So 

which is about 0(1000) times smaller than the value Yi ~ 2.5 x 10“^° which 
is required for the sufficient baryon number. Here, rii and Sq are the lepton 
number density and the entropy density today, and for simplicity, we neglect 
the SUSY contribution, which will be discussed later. 

An important observation for leptogenesis in this scenario is that under 
hxed Yukawa couplings, K oc 1/Ms^ and e oc Ms^. Therefore, larger Ms^ 
results in larger baryon number. This observation is critical because in this 
scenario, the mass of Si tends to be larger than expected by the symmetry. 
There are two essential points in this scenario. One of them is that it has 
a plenty of terms which give mass to Si. Each term gives the same order 
of mass to Si as expected by the symmetry, and the real mass can increase 
because of the large number of mass terms. The other point is that the 
predictions for the quark and lepton masses and mixings does not change so 
much even if the mass of Si becomes larger than expected by the symmetry. 
This is because the number of RH neutrino flavors becomes larger than three 
in Eg unihcation. (In S'O(IO) unihcation, it is not avoidable to change the 
predictions on neutrino sector if one of the RH neutrino masses is taken to 
be larger than expected by the symmetry.) 

The question is how large enhancement of the mass is needed to obtain 
the sufficiently large baryon number. It is the main subject in this paper to 
answer this question. 

In section 2, we briefly review the Eq GUT with anomalous U(1 )a gauge 
symmetry. And in section 3, we discuss the enhancement of the RH neutrino 
masses in this scenario. In section 4, we would like to answer the above ques¬ 
tion. First, we calculate the enhancement factor required to obtain the suffi¬ 
cient baryon number with simple non-SUSY Boltzmann equations. Second, 
we discuss the effect of lepton flavors. Third, the SUSY effect is considered. 
Finally, we calculate the enhancement factor, including both effects of lepton 
flavors and of SUSY. We show that only about 0(10) enhancement of the 
mass of Si is sufficient for the observed baryon number. 
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Table 1: Field contents of matters and GUT Higgs in a typical Eq X U{1U 
GUT[5] and the charge assignment nnder Eq x U{1)a- Here,Tj {i = 1,2,3) 
are three generation quarks and leptons, the VEVs of H and H break Eq 
into 50(10), the VEVs of C and C break 50(10) into SU{5), and the VEV 
of A breaks 50(5) into the standard model gauge group. The MSSM Higgs 
are included in H and C. 


2 Eq unification with anomalous U{1)a gauge 
symmetry 

We briefly review the EqxU{1)a GUT in this section [3]. The typical quantum 
numbers of helds in Eq x U{1)a GUT are shown in Table 1. An interesting 
structure in Eq unihcation is that three of six 5 of SU{5) in three matter 
helds \kj(27) become superheavy through the Yukawa interactions 

( 6 ) 

after developing the VEVs (HH) ~ and (OO) ~ which break 

Eq into 50(10) and 50(10) into 50(5), respectively. Here, the components 
of Yukawa matrices and are hxed by the total anomalous U{1)a 
charges of the corresponding terms TjTjTT and T^TjO, respectively. Since 
the Yukawa couplings for T 3 are larger than those for T 2 and Ti because 
V's ^ 53 and 83 become superheavy, and therefore, three light modes 

5 come from the and \k 2 . This structure naturally explains why 10s of 
SU{5) induce stronger hierarchy than 5s of 5f/(5), which is important to 
obtain realistic hierarchies of quark and lepton masses and mixings. 

The Eq X U{1)a GUT in Tabled] predicts the six RH neutrino masses 
(a = 1, 2, • • • , 6 ) and the Dirac neutrino Yukawa couplings Yai {i = 1, 2, 3) 
as in Table |2] except 0(1) coefficients. In the followings, we briefly review 
the derivation of these predictions from the model. See Ref. for the 
detail. The masses of the RH neutrinos can be obtained through the higher 
dimensional interactions 

(X, Y = H,C), (7) 

after developing the VEVs (H) ~ and (O) ~ A“ 5 F+'=). (These VEVs 

are determined by the VEV relations for the GUT singlet operators HH and 
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Table 2: GUT scale A^, Majorana 
1,2, and each component of nei 

GUT model with A = 0.22. 


Parameter 

value 

Ag 

2.000 X lO^^GeV 

Ml = A^^Ag 

5.656 X lO^GeV 

M 2 = A^^Ag 

2.571 X 10®GeV 

M 3 = A^^Ag 

1.169 X lO^GeV 

M 4 = A^oAg 

5.312 X lO^GeV 

Ms = A^Ag 

4.989 X 10^1 GeV 

Me = A^Ag 

2.268 X 10^2 GeV 

Yn = A6-5 

5.318 X 10-5 

Yi 2 = A 60 

1.134 X 10-^ 

Yi3 = A5-5 

2.417 X 10-^ 

Y 21 = A 60 

1.134 X 10-^ 

Y 22 = A^-s 

2.417 X 10-^ 

CO 

>- 

b 

5.154 X 10-^ 

Y 31 = A^'S 

2.417 X 10-^ 

Y 32 = A^-o 

5.154 X 10-^ 

>33 = A^-s 

1.099 X 10-3 

Y 41 = A^-o 

5.154 X 10-^ 

Y 42 = A^'S 

1.099 X 10-3 

^3 = A"-o 

2.343 X 10-3 

Y 51 = A3'5 

4.994 X 10-3 

^2 = A30 

1.065 X 10-2 

Ui3 = A2'5 

2.270 X 10-2 

Yei = A3-0 

1.065 X 10-2 

Y62 = A2-5 

2.270 X 10-2 

>63 = A20 

4.840 X 10-2 


masses of RH neutrinos (a = 
trino Yukawa U, in the Eq X C/(l).4 


comment 

GUT scale and the cutoff scale 

1st RH neutrino mass 

2nd RH neutrino mass 

3rd RH neutrino mass 

4th RH neutrino mass 

5th RH neutrino mass 

6th RH neutrino mass 

11 component of Y^, 

12 component of 

13 component of Y^, 

21 component of Y^, 

22 component of Y^ 

23 component of Y^y 

31 component of Y,y 

32 component of Y^y 

33 component of 

41 component of Y^, 

42 component of Y^, 

43 component of Y^, 

51 component of Y^, 

52 component of Y^, 

53 component of 

61 component of Y,y 

62 component of Y^y 

63 component of 



CC and the D-flatness conditions.) For example, the mass of Si 
becomes ;\2hi+2fe-(/j+/j)^ ^ ~ A^^A). It is convenient to 

dehne the effective U{1)a charges for any helds d' as 

^ = "0+^cy(d^) + (8) 


where cv'(tE') and cv/(\E') are the U{l)v and U{l)v' charges of respectively. 
The coefficients in the above eqnation ([8]) are determined so that the rela¬ 
tions {H) ~ A“^, {H) ~ A“^, (C) ~ A”'^, and {C) ~ A“^ are satished. It 
is obvious that the relations o and ([2]) do not change when the effective 
U{1)a charges are introduced because the U{l)v and U{l)v' charges of the 
Eq invariant terms are vanishing. Although special relations between 0(1) 
coefficients due to the Eq symmetry (or other original symmetries which 
are broken in the effective theory) cannot be seen explicitly in the effec¬ 
tive model, the effective U{1)a charges are useful to estimate the couplings 
of any terms allowed by the original symmetry. For example, the Dirac 
Yukawa couplings are easily estimated by these effective U{1)a charges as 
(^\h+n%.+hu^ ^ii+sj+hu-^_ The mass matrices of the RH neutrinos are also cal¬ 


culated as 


X^Ri+'^Rj \^Ri + ^3 

yi+n%j \Si + Sj 


The RH neutrino masses in Table H] 


can be obtained by diagonalizing the 6x6 RH neutrino mass matrix. In 
Table [2] we change the ordering of the RH neutrinos’ generation number a 
so that smaller number RH neutrino has smaller mass. The Dirac neutrino 
Yukawa couplings in Table [2] use this new index a. Even effective higher di¬ 
mensional interactions which give the light neutrino masses can be estimated 
as 


A 




UiH, 


(9) 


which are also derived from the RH neutrino mass matrix and the Dirac 
neutrino Yukawa matrix by the seesaw mechanism. 

One of the most interesting features in the anomalous U{1)a models is 
that the higher dimensional interactions give the same contributions to in¬ 
teractions as the lower dimensional interactions. For example, the coeffi¬ 
cients of Yukawa interactions TjTji/ are determined by their total U{1)a 
charge as except 0(1) coefficient. The higher dimensional interac¬ 

tions TjATjiif, whose coefficients are also determined by the total charge as 
also contribute to the Yukawa interactions TjTjFf after devel¬ 
oping the VEV (A) ~ A““ which breaks SU{5) into the SM gauge group. 
The coefficients from the higher dimensional interactions are estimated as 
\ilji+ilj 3 +h+a ~ which is nothing but the coefficients of the orig¬ 

inal Yukawa interactions except 0(1) coefficients. Therefore, the unrealistic 
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GUT relations of Yukawa couplings, for example, can be naturally 

avoided in the anomalous U{1)a GUT models because the higher dimen¬ 
sional interactions with the adjoint Higgs A have different contributions to 
the down-type Yukawa couplings from the charged lepton Yukawa couplings 
after developing the VEV of A. 

3 Possible enhancement for the right-handed 
neutrino masses 

It is plausible to enhance a coefficient of an interaction if there are a lot 
of higher dimensional interactions which contribute to the coefficient by the 
same order after developing the VEVs of the negatively charged operators. 
Roughly, if there are N higher dimensional interactions which give the same 
contribution to an interaction, the enhancement factor can be expected to be 
y/N according to the random walk theory. Since we have introduced several 
negatively charged singlets as well as the GUT Higgs helds, the number N 
can be large if the total U{1)a charge of an interaction is large. For example, 
in a simplihed model in which all negatively charged helds ©» (z = 1, 2, • • • , n) 
have the U{1)a charges 6i = —1, the number of the independent interactions 
with total U{1)a charge c is given by Nn{c) = • This number Nn{c) 

becomes easily large when c and n are large. For example, we obtain that 
^ 5 ( 5 ) = 126, iV5(10) = 1001, Yio(lO) = 92378, ■ ■ ■. In this section, we will 
show it is plausible that the 1st, 2nd and 3rd smallest RH neutrino masses 
are enhanced and this enhancement does not change the physical predictions 
for the light neutrino sector so much. 

The interactions which contribute to the masses of the RH neutrinos Si 
and (z = 1, 2, 3) are ^A^iHH and TjTjC'G, respectively. The total U{1)a 
charges of these interactions are (11, 9, 5) for Si and (9, 7, 3) for while 
the masses expected by the symmetry are (A^^, A’^) and (A^^, A^°, A®), 

respectively. This means that the enhancement factors r|s^ and for their 
masses are expected to be the largest for the lightest RH neutrino Si, the 
second largest for the second and the third lightest neutrinos and 5*2. 

In this paper, we do not count the total number of the independent in¬ 
teractions which give the mass term of these RH neutrinos in the explicit Eq 
GUT model in Table [H However, we discuss what happens when some of 
the RH neutrinos have larger masses than those expected by the symmetry. 
It is an important observation that each RH neutrino gives the same order 
of the contribution to all components of the light neutrinos’ mass matrix 
My = Y^^M~^Yy^{Huy if its mass is nothing but the value expected by 



the symmetry. Therefore, if one of the enhancement factors r|s^ and is 
aronnd one, all components of M^, becomes the valnes expected by the sym¬ 
metry, and so are all components of the diagonalizing matrix. In order to 
obtain three eigenvalues expected by the symmetry, three of the six enhance¬ 
ment factors must be around one. Then all predictions on the light neutrino 
sector become the same order as the predictions without the enhancement 
factors. Since the lightest neutrino mass has been fixed only its upper limit 
by experiments, the prediction for it can be different from the predicted 
value without any enhancement factor. Therefore, it is sufficient that two 
RH neutrinos have their masses which are determined by the symmetry for 
consistency with the present constraints obtained by neutrino experiments. 

It looks not to be fair that we consider these enhancement effects only for 
the RH neutrino masses, although the mass terms have much larger U{1)a 
charges than the other terms like Yukawa terms. We should change the U{1)a 
charge assignment in Table [U when such enhancement effects are taken into 
account. This subject is beyond the scope of this paper. Here we should 
emphasize that even after changing these U{1)a charges, the mass terms of 
Si, S 2 , and have still much larger U{1)a charges, and therefore, some 
enhancements for their masses are expected. 

The next important question is how large enhancement factor is needed 
for sufficient leptogenesis in this Eq GUT model. In the next section, we try 
to answer this question. 


4 Leptogenesis in the Eq x U{1)a model 

In the thermal leptogenesis scenario, thermally produced RH neutrinos go 
out of equilibrium as temperature decreases to their mass scale, and their CP 
asymmetric decays produce lepton asymmetry [12] . The lepton asymmetry is 
converted to the baryon asymmetry via the nonperturbative B + L violating 
sphaleron processes HI. 

In this section, we calculate the thermally produced lepton number in 
the Eq X U{1)a model with the Dirac neutrino Yukawa couplings Y^i {a = 
1, 2, • • • ,6,i = 1,2,3) which are determined by the symmetry as in Table [2] 
and the masses Ma for the mass eigenstate of the RH neutrinos Na. Some 
of six Ma have enhancement factors rja larger than 1. What we would like 
to know by this calculation is how large enhancement factors are required to 
obtain sufficiently large lepton number. In the calculation, it is important to 
include supersymmetric contributions and the effects of lepton flavor in the 
hnal state of the decay process simultaneously. To show this statement, we 
calculate the sufficient enhancement factor in four cases: 
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Figure 1: r/i = Mi/M° dependence of |Fb-l| in each case. Horizontal band 
corresponds to the observed baryon asymmetry in SUSY cases. is the 
“bare” Majorana mass in the absence of the U{1)a interactions. We take the 
simplihed CP asymmetry e™ ffTHjl . e^usY = 2 x e™, and ^ 

2 X with the assumption Q^[(YW)6 i] = 3fJ[(YW)6i] for the calculation in 
each case, respectively. and i^susY written as ~ ST/r^i and 

/^susY ^ 51 / 771 , respectively. 


• non-SUSY + non flavor 

• non-SUSY + flavor 

• SUSY + non flavor 

• SUSY + flavor 

The result is shown in Fig. [TJ In a realistic situation of the Eq x U{1)a GUT 
model, i.e., in the case of SUSY-1-flavor, the sufficient lepton number can be 
obtained if the enhancement factor for the Ni mass is around 16. This means 
that Ml ~ 9 X 10® GeV. 

It is known that supersymmetric contribution is important when K is 
smaller than 1 because supersymmetric calculation makes K larger effec¬ 
tively. On the other hand, the lepton flavor effects are important when the 
decay parameter K = Fiv^(T = 0)/H(T = Mi) is larger than 1 because K 


10 









Figure 2: Tree and one-loop diagrams contributing to the CP asymmetric 
decay of lightst RH neutrino. In the model, hve RH neutrinos Np (/3 = 
2,6) contribute to the CP asymmetry. 


for the muon and the electron become smaller than K for the tan. Here T 
is temperature of the universe. Our calculation shows that it is important 
to inculude both contributions when K is larger than 1. This is because 
supersymmetric contribution is important for smaller K of the electron and 
the muon. 

4.1 non-SUSY + non-flavor 

In this subsection, we evaluate the lepton asymmetry in the non-SUSY-|-non- 
flavor case in Eq x U{1)a models after brief review. This simple calculation 
is important to understand the outline of the thermal leptogenesis with the 
RH neutrino masses and the Dirac Yukawa couplings in Table [2l 

In the model, since RH neutrinos are hierarchical in mass, the lepton 
asymmetry is generated by the CP asymmetric reactions of the lightest RH 
neutrino Ni. In the followings, we assume that the lightest RH neutrino is 
Ni while it has the largest enhancement factor. The lepton asymmetry is 
evaluated by a coupled set of evolution equations of the lightest RH neutrino 
Ni and the lepton asymmetry L\ 
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Figure 3: Evolutions of \Yb-l\) Ynd and for Mi = 1 xMj* (top panel) and 
Ml = 100 X (bottom panel) in the non-SUSY+non-flavor case. Here 
is the “bare” Majorana mass in the absence of the U{1)a interactions. Hori¬ 
zontal band corresponds to the observed baryon asymmetry. We take the sim- 
plihed CP asymmetry flT^ with the assumption 9^[(Y1Y)6 i] = 3f?[(Y'^Y)6i]. 


dYL 

dz 


sH{z = 1) { [2 Y;"'' 


1 Y 


L I ^SM 
■r ^Ni 


Yj 


Vi 




Id + 


Yl 

Yr i 


9 ^ 

^'ySt Y -.^en'ySs 




( 11 ) 


Here s is the entropy density, and H is the Hubble parameter. We use 
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a dimensionless variable 2 ; = Mi/T. Yx and are yield value and its 
equilibrium one of a species X, respectively, which are the number density 
normalized to the entropy density. 7 ^ = 7 d(Xi o IH), jss = 1 Ss{Nil ■(-)■ 
QsT^), and 'jst = IstiNiQ^ ■(-)■ ITr) are thermal averaged decay rate (inverse 
decay rate), s-channel, and t-channel scattering rate, respectively [22]. Here 
/, Qs, and Tr denote SU{2)r doublet lepton, third generation doublet quark, 
and singlet top quark, respectively. The CP asymmetry is dehned as 
eff = [r(Xi ^ IH) - r(iVi ^ IH) + r(Xi ^ IH% The 

hrst non-zero contribution to comes from interference between tree-level 
amplitude with the one-loop contributions (upper three diagrams in Fig. [2|), 
and it is calculated in a hierarchical limit in RH neutrino masses as = 

-(3/167r)E?^i (Mi/M^) [I5]. Note that, the E, 

GUT model has six RH neutrinos, and therefore, /3 = 2, 3, • • • , 6 . 

Key ingredients for the lepton asymmetry generation are the CP asym¬ 
metry and the decay parameter K = = 0)/H{T = Mi) [121 

which parametrizes the departure from the thermal equilibrium of RH neu¬ 
trinos at T = Ml. K is important because it is related with 7 ^) and the 
factor (1 — Yxi/Y^) in Eq. ffTTjl . The lepton asymmetry Yl is essentially 
determined by the above two parameters as 

Yl ~ 4“C(A'). (12) 

The behavior of the function C{K) is as follows. When K > 1, C{K) becomes 
a decreasing function of K. K > 1 means that the RH neutrinos are still 
in the thermal equilibrium at T = Mi, and therefore, the number density 
of Ni decreases rapidly when T < Mi. This reduces the produced lepton 
asymmetry. Obviously larger K results in lower decoupling temperature, 
smaller Yxi after the decoupling and smaller lepton asymmetry. When K < 
1, C{K) becomes a increasing function of K. K < 1 means that the RH 
neutrinos are out of the thermal equilibrium at T = Mi, and therefore, the 
number of thermally produced RH neutrinos becomes smaller for smaller K. 
This reduces the produced lepton asymmetry. Around iP ~ 1 , the function 
C{K) becomes maximal. Sufficient lepton asymmetry can be obtained when 
A ~ 1 and e™ ~ 10 “’^. 

Let us calculate the above two important parameters in the Eq XC/( 1 U 
GUT. First, we estimate and without any enhancement factor rja 
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as 


^SM 


^SM 


2 f 3 ^[{Y^Y)1] mA 

\ 167r Me) 

[Y^Y]uM,/87r ^ 


-8.85 X 10”® 


5.7 X 10^GeV\ 

) 


Ml 


5.7 X 107 GeV 



(14) 


where Mpi is the Planck scale and is the effective relativistic degrees 
of freedom, which is obtained as gf^ = 106.75 with the SM particle con¬ 
tents. For the estimation of e^, we have adopted two assnmptions. The 
first assnmption is that can be estimated by 3?. This 

assnmption is reasonable because we regard all Yukawa couplings as complex 
numbers. The second assumption is on the factor 2 in front of the parenthesis 
in Eq. flT^ . An important observation is that ^ jg 

not dependent on /3. Therefore, we can expect an enhancement factor after 
summation of the index /9, and we assume that the enhancement factor is 
two through all calculations in this paper. 

It is obvious that the lepton asymmetry with these parameters are too 
small to explain the observation. is too large and e|^ is too small. How¬ 
ever, as discussed in the previous section, the lightest RH neutrino mass can 
be expected to have an enhancement factor which can be much larger than 
one. Interestingly, when the lightest RH neutrino mass Mi becomes larger, 
the produced lepton asymmetry becomes larger because the GP asymme¬ 
try becomes larger and the decay parameter becomes smaller as 
seen in Eq. flTSl) . For example, if we take the enhancement factor is around 
37, the sufficient lepton asymmetry can be expected because ^ 

~ 3 X 10-7. 

Figure [3] shows the evolutions of the lepton asymmetry \Yb-l\, yield value 
of the RH neutrino Yiv^, and its equilibrium one Y^^ for Mi = 1 x Mj* 
(top panel) and Mi = 100 x Mj* (bottom panel). Here M° represents the 
“bare” Majorana mass, that is the physical mass of the lightest RH neutrino 
without any enhancement factor. The lepton asymmetry for Mi = 1 x M°, 
\Yl\ ~ 10“^^, is too small to account for the observed baryon asymmetry. In 
non-SUSY cases, the required lepton asymmetry is 2.285 x lO-^^ < Yb-l ^ 
2.685 X 10-7° with the conversion rate of the lepton asymmetry to the baryon 
asymmetry Yb = {28/79)Yb-l [E] and the observed baryon number 8.097 x 
10-77 <Yb< 9.518 X 10-77 [T7]. For Mi = 100 x M°, the lepton asymmetry 
is drastically enhanced. The enhancement of the lightest RH neutrino mass 
makes the GP asymmetry larger [see Eq. (IT^ ] and reduces the factor. 
Indeed, in the bottom panel in Fig. [31 we find the larger deviation from 
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thermal equilibrium compared with the top panel. The combination of these 
effects leads the enhancement of lepton asymmetry. 

Dotted line in Fig. [T]shows the Mi dependence of the lepton asymmetry in 
the non-SUSY+non-flavor case. Since numerically C{K) oc K for iF 1, the 
lepton asymmetry becomes asymptotically a constant for the enhancement 
factor larger than 40. In this case, the physical mass of the RH neutrino 
is required to be Mi = (35 — 39) x to account for the observed baryon 
number. 


4.2 non-SUSY + flavor 


It is important for the evaluation of lepton asymmetry in the Eq X U(1U 
model to separately involve each lepton flavor channel of the CP asymmetric 
decays. The reasons are as follows. The Eq model possesses the features: 
(i) the evolution of lightest RH neutrino is in the strong washout regime, (i 
i) all of asymmetry productions of each lepton flavor by the CP asymmetric 
decay are sizable. These features can give rise to 0{1) corrections to the 
final lepton asymmetry with respect to the case where the flavor effects are 
ignored ra. This is because that the evolutions of the lepton asymmetries 
of each lepton flavor are in the regime of washout with different magnitudes. 
In this section we briefly review the leptogenesis with the flavor effects. 

The total lepton asymmetry is given by the sum of the asymmetry of each 
lepton flavor, Yb-l = pAe + + Pa.^) where Aj = B/?> — L^. The subscript 

i represents the lepton flavor. The evolutions of each asymmetry is described 
by the flavor dependent Boltzmann equations. 


dY, 


Ai 


dz sH{z = 1) 

'Pv, 


Yj 


Ni 




1 4“7« + /y^ i(4 + c«)7D 




CH 


+ ( ^ - 1) I CljlSs + -^ISt ) + (2^^ + cf) (ys, + ^ 


ISs 


yeg 


(15) 


The coefficient is the ffavor projection, = Yi 4 Y'i’i/(P"P^^)ii- The ffa- 
vor dependent CP asymmetry is defined as eff' = [P(iVi —)■ liH) — T{Ni —)■ 
liH^)]/[T{Ni —)■ liH) + P(iVi liW)], and is calculated in the hierarchical 


mass limit as = - {l/8n{YY^)ii) Yp^Y* (3/2) (Mi/M^) (YRt) 

+(m2/m 2) (yyt) 


1/3 


'/31 


T5] . We follow the considerations for deriving Eq. m, 


and obtain the simplihed CP asymmetry in the non-SUSY+flavor case as 


^SM _ 9 
^ 




-1 


V87r(RRt) 
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( 16 ) 
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Figure 4: Same as Fig. |3l but in the non-SUSY+flavor case. The evolutions 
of electron number L^, muon number and tan number are also plotted. 
Since the evolution of RH neutrino is same as the non-SUSY+non-flavor case, 
we cut the part. 

The coefficient O'" {C^) in Eq. f[T5|l is introduced as the conversion factor 
between the asymmetry normalized to the equilibrium number density for 
(H) and the yield value of each lepton number normalized to equilibrium lep¬ 
ton density as {rii. - njJ (Fh , and {hh - Uh) /n^ = 

— entries are model-independent, which are deter- 
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mined by constraints among the chemical potentials enforced by the equilib¬ 
rium reactions in the temperature T ~ Mi where the asymmetries are gener¬ 
ated. The region of RH neutrino mass we consider is 1 < Mi/M° < 100, and 
the relevant temperature of the leptogenesis in the model is in 10® GeV < 
T < 10^^ GeV. In this range, SM gauge interactions, third and second gener¬ 
ation Yukawa interactions are in equilibrium, and the equilibrium conditions 
lead to the following and [20] , 



906 

-120 

-120\ 

1 


-75 

688 

-28 , 

C^-gSgT 52 52). 

(17) 

-75 

-28 

688 / 



Figure m shows the evolutions of total lepton asymmetry \Yb-l\^ and of 
the asymmetry of each lepton flavor | YaJ (* = e, /i, r) for Mi = 1 x (top 
panel) and for Mi = 100 x (bottom panel). For the interpretation of the 
result, we need to see both the magnitude of the washout and the production 
efficiency of each lepton asymmetry. We rearrange the factor ffTTl) to 
involve the flavor dependence as Kf^ = Each Kf^ is obtained as 

follows. 


K: 


SM 


= 0) ^ / 5.7xlO^GeV 
H{T^Mi) ' f Ml 


( 18 ) 


SM = 0) , / 5.7 X 10^GeV 

^ H{T = Mi) ■ V Ml 

'SM ^ = 0) ^ lyq f 5-7 x 10^ GeV 

HiT = Mi) - \ Ml ^ 


(19) 

( 20 ) 


The Kf^ is a measure of magnitude of the washout of each lepton asym¬ 
metry, that is the same with the relation between the factor and the 
washout of total lepton asymmetry. On the other hand, the ratio of asym¬ 
metry productions of each lepton flavor by the GP asymmetric decay is equal 
to the ratio of each GP asymmetry, and is obtained as follows 


Le production ^ ^ [> 61 ^ 1 * 1 ] ^ ^ 

Lr production 3 [Vishlg] ‘ ’ 

Lf, production ^ ^ [> 62 ^ 1 * 2 ] ^ ^ 220 
Lr production [^ 63 ^" 13 ] 

Here we assumed O' [Vei^n] = 3^^ For Mi = 1 x M^, and are 

(P(l), but Kr — 29. Each lepton asymmetry is generated by the GP asym¬ 
metric decays, and then a large part of the tau’s is washed out by the inverse 
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decays and so on, while the electron’s and the mnon’s soon deconple from 
the eqnilibrinm and snrvive. Thns, nonetheless the prodnction efficiencies of 
the electron and the mnon nnmber are lower than the tan’s, they yields a 
large part of the lepton asymmetry. On the other hand, for Mi = 100 x M°, 
Kr — 0(10“^), and ~ 0(10“^). These Ki factors indicate that each 

lepton nnmber generated by the CP asymmetric decay snrvives withont be¬ 
ing strongly washed ont. Thus a large part of lepton asymmetry is governed 
by the tan’s. 

The chain line in Fig. [T] shows the Mi dependence of the total lepton 
asymmetry in the non-SUSY-|-flavor case. In the non-SUSY-|-flavor case, 
the physical mass of the RH neutrino is required to be 30 < Mi/M^ < 37 
to account for the observed baryon number. For Mi/M° < 40, in the non- 
SUSY-I-non-flavor case, the evolution of the lepton asymmetry is in the strong 
washout regime. While, in the non-SUSY-|-flavor case, the muon and the elec¬ 
tron asymmetries are not hardly washed out, and yield sizable contribution 
to total lepton asymmetry. The lepton asymmetry with the flavor effects is 
therefore larger than that in the case where the effects are ignored. On the 
other hand, for Mi/M° > 40, the evolutions of the asymmetries of all lepton 
flavor are in the weak washout regime, and hence total lepton asymmetry is 
determined by only the asymmetry production by the CP asymmetric decay. 
Due to the additional washout contributions, in the parameter region, the 
hnal lepton asymmetry can be smaller than that the case without the flavor 
effects [ISl 120] • 

4.3 SUSY + non-flavor 

The SUSY extension of the leptogenesis gives an enhancement for the lepton 
asymmetry which is roughly estimated as 

yjusY I' ^ (strong washout) 

Y™ ^ \ 2^/2 (weak washout) 

in Ref. [ 21 ] • In this section, we briefly review the corrections to interpret the 
numerical results in the context of the Eq x U{1)a model. 

We have two important points which increase the lepton asymmetry. The 
additional decay channels correct the dehnition of the CP asymmetry as 
e^usY _ ^y{Ni lW) + T{m ^ /M)-F(Yi ^ PM)]/F^usy_ 

Similarly, for the RH sneutrino, = [r(iVi IH) — F(iVi —)■ Pifl) -|- 

F(iVi ^ IH) - F(iVi ^ rM)]/F^usY_ YLeie F^^sy jg total width of the 
RH sneutrino. The CP asymmetry receives the contributions of not only 
the RH neutrinos also of its scalar partner, and is obtained as = 
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— (3/87r) (^ 1 /^/ 3 ) the hierarchical limit of 

RH neutrino masses na. The CP asymmetry of the RH sneutrino is equal 
to that of the RH neutrino in the hierarchical mass limit. Repeating the 
consideration for deriving Eq. ca, we obtain the simplihed CP asymmetries 
in the SUSY+non-flavor case as follows, 


^SUSY _ ^SUSY 


3 mY^Y)l] MA 
Stt M,) 


2 X e' 


SM 
iVi • 


(23) 


These effects make the lepton asymmetry four times larger. 

However, the effective relativistic degrees of freedom = 228.75 is 

about twice of which reduces the lepton asymmetry to entropy ratio by 

1/2. 

The most important one is the correction of K factor. In the context of 
SUSY, there exists additional decay channels of the RH neutrino, and the 
total width is obtained as P^®"^ = r(Yi —)■ IH) + P(Yi —)■ IH^) + P(Yi —)■ 
IH) + P(iVi ^ l*H) = [YtY]nMi/47r = 2 x P^^. I and H represent SUSY 
partners of the SU{2)l lepton and the Higgs doublet, respectively. The K 
factor is calculated as follows, 

,(-susY _ = 0) _ lKty|iiAfi/4.r ^ / 5.7x lO^GeV A 

H{T = Ml) 1.66(^susY)i/2jvf2/j\^^j j ■ 

(24) 


Roughly, the factor iySUSY jg times larger than that in the SM. This effect 
reduces the lepton asymmetry in strong washout regime and enhances it in 
weak washout regime. 

Finally we note the conversion rate from the lepton asymmetry to the 
baryon asymmetry. In the context of SUSY, the additional equilibrium reac¬ 
tions at the temperature T ~ Mi alter the constraints among the chemical 
potentials. The alteration leads the conversion rate as Yb = {9>/2?))Yb-l [21]- 
Consequently the required lepton asymmetry in SUSY cases is 2.328x10“^° < 
\Yb-l\ < 2.736 X 10-10. 

With all these effects, the result in Eq. fl22]) is obtained. 

Dashed line in Fig. [1] shows the Mi dependence of the lepton asymmetry 
in the SUSY-1-non-flavor case. The lepton asymmetry is given by a sum of 
partial asymmetries from the CP asymmetric decays of RH neutrino and its 
scalar partner. It is evaluated by a coupled set of evolution equations of the 
RH neutrino, its scalar partner, and the partial asymmetries [22]. Due to too 
strong washout, for Mi < 30 xM)*, nonetheless the additional CP asymmetric 
decays, the lepton asymmetry is close to that in the non-SUSY-1-non-flavor 
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Figure 5: rji = Mi/M° dependence of |Fb-l| and IIaJ in the SUSY+flavor 
case. Horizontal band corresponds to the observed baryon asymmetry. We 
take the simplified CP asymmetry = 2 x with the assumption 

Q^[(yt’F)6i] = 3fJ[(yty)6i], where e™ is given by ffTHjl . can be written 

as i^SUSY _ 


case. While, for Mi > 30xM°, because of both K ~ 0(1) and the additional 
CP asymmetric decays, larger lepton asymmetry is generated than those in 
the cases of non-SUSY+non-flavor and non-SUSY+flavor. 

4.4 SUSY + flavor 

We are now in a position to discuss the lepton asymmetry in the SUSY+flavor 
case, which involves the SUSY particles contributions with the flavor effects 
[see Sec. K2\. This is the realistic situation in the Eq x U{1)a model. In¬ 
terestingly, even in the strong washout regime which is dehned by > 1, 
the effect of SUSY becomes sizable if the flavor effects are included. 

The sum of total lepton and slepton asymmetries, and , converts 
to the baryon asymmetry, and which are given by the sum of the asymmetry 
of each lepton and slepton flavor: Yb-l = Y^^^ + Y^^'^ = (Yab + + Y^^) + 

(Y^^ + Y^^ + Y^^). We take the SUSY spectrum to be 0(1 TeV) in the 
EqxU (1)a model. Then, throughout the temperature region we consider, the 
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Figure 6: Evolutions of \Yb-l\ and |FaJ for Mi = 16 x in the 
SUSY+flavor case. Horizontal band corresponds to the observed baryon 
asymmetry. We take the simplihed CP asymmetry = 2 x with the 

assumption ^>[(1^11^)61] = 3fJ[(y'^Y)6i], where e™ is given by ffTHjl . 


equality of chemical potentials of a SM particle and its superpartner, which is 
refferd as superequilibration |2S], is maintained in the presence of equilibrium 
supergauge (or Yukawa) reactions. In the superequilibrium regime, since 
the equality of asymmetries of each lepton and its scalar partner is also 
maintained, Yb-l = 2 x (Ya^ + Thus the baryon asymmetry 

in the case is evaluated by a coupled set of evolution equations of the RH 
neutrino, its scalar partner, and asymmetry of each lepton flavor. The flavor 
dependent Boltzmann equations are shown in Ref. [23], and relevant cross 
sections are given in Ref. |22j . 

Figure [5] shows the Mi dependence of the lepton asymmetry |Yb-l| and 
the partial asymmetry of each lepton flavor |YaJ- From the Fig. [5l the en¬ 
hancement factor 16 < Mi/Mf < 17 can yield the observed baryon number. 
The Ml dependence of the asymmetries are described by the washout effects 
and production efficiencies of the asymmetries, that is basically the same 
as in the non-SUSY-|-flavor case. The magnitude of washout of each lepton 
flavor is parametrized by which is SUSY extension of Kf^ [Eqs. fflSj) 
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Figure 7: rji = Mi/Mf dependence of the ratio of flavor con¬ 
tributions in each case, (|1 aJ)susy / (I^aJ)sm> (I^b-^DLsy / 
(l>^B-n|)rus-Y/(l>^B-n)rj- ^ Here, for example, |Fb-l|sm lepton asym¬ 

metry in the case non-SUSY-|-flavor, and IFb-lIsu'sy^ is the lepton asymmetry 
in the case SUSY-1-non-flavor, etc. i^susY written as iySUSY ^ 5 I/ 771 . 


- m]-- 


rSVSY 


= 0) ^ / 5.7xlO^GeV 

H{T^Mi) ' 1 Ml 


(25) 


SUSY = 0) ,, 1 , 1 , / 5.7 X 10^ GeV 

^ H{T = Mi) ■ V Ml 

'SUSY ^ = 0) ^ SQ / 5.7 X 10^ GeV 

H{T = Mi) - \ Ml ^ 


(26) 

(27) 


As in the non-flavor case, these K factors become \/2 times larger than in 
the non-SUSY case. These corrections make the washout of each asymmetry 
stronger compared with non-SUSY case and the SUSY effects become weak 
especially in the strong washout regime. However, even if > 1, some 
of the can be smaller than one, and therefore, the washout effect for 
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the flavor i is negligible. Then the supersymmetric contribution become siz¬ 
able. Consequently, the lepton asymmetry generation is sufficiently boosted 
compared with the case of the non-flavor especially when > 1. 

Figure [6] shows the evolutions of total lepton asymmetry |Fb-l| and par¬ 
tial asymmetries of each lepton flavor IFaJ for Mi/M° = 16. To understand 
the importance of the flavor effects in SUSY calculation, we make Figure [7] in 
which the ratios of SUSY lepton asymmetry to non-SUSY lepton asymmetry 
are plotted. First of all, the figure shows that SUSY enhancement factor 
is larger in the weak washout regime than in the strong washout regime as 
explained in the previous subsection. Next, the figure shows that the SUSY 
enhancement factor with the flavor effect is larger than without the flavor 
effect. Especially, it is important that even in the strong washout regime, 
the SUSY enhancement factor |YB_i|gugY/|ks-LlsM become sizable due to 
the enhancements of the muon and electron asymmetry, because |Yb-l|sm 
has already been fairly larger than |Yb-l|sm~^ strong washout regime. 

We could confirm the successful baryon asymmetry in the SUSY-|-flavor 
case, namely, in a realistic situation of the Eq x U{1)a_ GUT model. It is 
sufficient to take the lightest RH neutrino mass Mi ~ 16 X ~ 9 X 10® GeV 
for the observed baryon asymmetry. It is important for this calculation that 
all components of neutrino Yukawa matrix are determined by the symmety 
in the Eq x U{1)a GUT and we can integrate the flavor effects on the lepton 
asymmetry. 


5 Summary and Discussion 

We have investigated the thermal leptogenesis in the Eq x U{1)a GUT model 
in which realistic quark and lepton masses and mixings are obtained and the 
doublet-triplet splitting problem is solved with natural assumption that all 
interactions including higher dimensional interactions are introduced with 
0(1) coefficients. Each of three fundamental representations 27 includes two 
SM singlet fields, *S'(1’) and Y)j(l), and these singlet fields play a role of 
RH neutrinos (a = 1, 2,..., 6). One of the aim of this work is to show a 
sufficient lepton asymmetry is generated by the GP asymmetric decays of the 
lightest RH neutrino. In the model, Majorana masses of the RH neutrinos 
Ma and the neutrino Yukawa couplings Y^i are determined by the U{1)a 
symmetry. So we can calculate the lepton asymmetry, but unfortunately the 
naive calculation results in too small abundance of the lepton asymmetry. 
Actually, the lightest RH neutrino mass is around 6 x 10^ GeV, which is 
smaller than the Ibarra’s lower bound 10®“® GeV. Moreover, the factor K 
and the GP asymmetry are evaluated as iP ~ 40 and ~ 5 x 10“®. 
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Therefore, the lepton asymmetry is washed out strongly in this scenario, and 
even with /T ~ 1, the is too small to obtain the sufficient number of 
lepton asymmetry. 

We have shown that a key ingredient for successful leptogenesis is the 
enhancement of RH neutrino masses. The model can include a large number 
of higher dimensional interactions which, and these interaction terms yield 
additional Majorana masses after developing the VEVs of negatively U{1)a 
charged helds. The enhancements of the RH neutrino masses enhance the 
CP asymmetry e oc Mi and make the decay parameter K oc 1/Mi smaller 
to be most efficient value iP ~ 1. How large enhancement factor is required 
for the sufficient leptogenesis? To answer this question, we have calculated 
the lepton asymmetry including the effects of SUSY and flavor in the hnal 
state of the CP asymmetric decay. The result is that the enhancement factor 
16 —17 is sufficient for the successful leptogenesis. About 300 mass terms are 
sufficient to obtain this enhancement factor, and this number looks not to be 
difficult to be obtained even in the Eq GUT model. It is important that such 
enhancement of the lightest RH neutrino mass does not change the neutrino 
physics at the low energy scale. This is because the Eq x U{1)a GUT has six 
RH neutrinos which induces the same order of the amplitude of all elements 
of the light neutrino mass matrix. 

We have calculated the lepton asymmetry in the Eq x U{1)a model in 
following four cases: (i) non-SUSY+non-flavor (ii) non-SUSY+flavor (iii) 
SUSY+non-flavor (iv) SUSY+flavor. These calculation has shown that both 
the effects of lepton flavor and SUSY are important. It is known that in the 
strong washout regime lepton flavor effect becomes sizable, though SUSY 
contribution is not so large. We have shown that SUSY contribution be¬ 
comes important even in the strong washout regime if lepton flavor effect is 
included. The essential point is that even in the strong regime > 1, 
the washout effects of the muon and/or the electron can become weak, and 
therefore these lepton number abundances become sizable. 

Of course, the obtained result for the enhancement factor 16 — 17 (Mi ~ 
9 X 10® GeV) for the sufficient leptogenesis is dependent on the various pa¬ 
rameters and even on the 0(1) parameters. For example, we have hxed the 
coefficient of in Eq. f[T51) as two in our calculation. Since the hnal lep¬ 
ton asymmetry is proportional to this parameter, the dependence can 
be read from the Fig. [H When the coefficient is one, the required enhance¬ 
ment factor becomes around 25 (Mi ~ 1.4 x 10® GeV), and therefore, 0(600) 
mass terms are needed. When the coefficient is four, it becomes around 10 
(Ml ~ 6 X 10® GeV), which is required 0(100) mass terms. Therefore, we will 
not predict the mass of the lightest RH neutrino for sufficient leptogenesis, 
because it depends on the various parameters. An important thing is that 
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the Eq X U{1)a gut can explain the baryon asymmetry in the nniverse. 

It is not plansible to prodnce the snfiicient lepton nnmber in the S'O(IO) x 
U{1)a GUT[7] by the enhancement of the lightest RH nentrino. Since the 
nnmber of the RH nentrinos is three in S'O(IO) model, the other two nentrinos 
mnst have the masses expected by the symmetry. However, the difference 
between the U{1)a charges of the lightest RH nentrino and the second lightest 
RH nentrino mass terms is just two, and therefore, it is not reasonable to 
expect that the lightest RH neutrino has 0(10) times larger mass than the 
second lightest RH neutrino. However, since the Eq x U{1)a has six RH 
neutrinos and the difference between the U{1)a charges of the lightest RH 
neutrino and the forth lightest RH neutrino mass terms is six, it is plausible 
that the lightest RH neutrino has 0(100) times larger number of mass terms 
than the forth lightest RH neutrino. The observed baryon asymmetry in our 
universe may be an indirect signature of Eq GUT. 
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